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Abstract —Millimeter wave (mmWave) systems will likely em¬ 
ploy large antenna arrays at both the transmitters and receivers. 
A natural application of antenna arrays is simultaneous transmis¬ 
sion to multiple users, which requires multi-user precoding at the 
transmitter. Hardware constraints, however, make it difficult to 
apply conventional lower frequency MIMO precoding techniques 
at mm Wave. This paper proposes and analyzes a low complexity 
hybrid analog/digital precoding algorithm for downlink multi¬ 
user mm Wave systems. Hybrid precoding involves a combina¬ 
tion of analog and digital processing that is motivated by the 
requirement to reduce the power consumption of the complete 
radio frequency and mixed signal hardware. The proposed 
algorithm configures hybrid precoders at the transmitter and 
analog combiners at multiple receivers with a small training and 
feedback overhead. For this algorithm, we derive a lower bound 
on the achievable rate for the case of single-path channels, show 
its asymptotic optimality at large numbers of antennas, and make 
useful insights for more general cases. Simulation results show 
that the proposed algorithm offers higher sum rates compared 
with analog-only beamforming, and approaches the performance 
of the unconstrained digital precoding solutions. 

I. Introduction 

The large bandwidths in the mmWave spectrum make 
mmWave communication desirable for wireless local area 
networking and a candidate for future cellular systems Q- 
Achieving high quality communication links in mmWave 
systems requires employing large antenna arrays at both the 
access point or base station (BS) and the mobile stations 
(MS’s) Q, 0. The transmit antenna array can be used to 
send data to multiple users using multi-user MIMO pre¬ 
coding principles. In conventional lower frequency systems, 
this precoding is usually performed in the digital baseband 
to have a better control over the entries of the precoding 
matrix. Unfortunately, the high cost and power consumption 
of mixed signal components make fully digital baseband pre¬ 
coding difficult at mmWave systems 0. Further, constructing 
precoding matrices is usually based on complete channel state 
information, which is hard to acquire in mmWave systems 
0. Therefore, mm Wave-specific multi-user MIMO precoding 
algorithms are required. 

In single-user mmWave systems, analog beamforming, 
which controls the phase of the signal transmitted at each 
antenna via a network of analog phase shifters and is im¬ 
plemented in the radio frequency (RF) domain, is the defacto 
approach for implementing beamforming 0. The phase 
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shifters might be digitally controlled with only quantized phase 
values. In 0, adaptive beamforming algorithms and multi¬ 
resolution codebooks were developed by which the transmitter 
and receiver jointly design their analog beamforming vectors. 
In 0, beamspace MIMO was introduced for large antenna 
systems in which discrete Fourier transform beamforming vec¬ 
tors are used to asymptotically maximize the received signal 
power. The RF hardware constraints such as the availability 
of only quantized and constant modulus phase shifters, how¬ 
ever, limit the ability to make sophisticated processing using 
analog-only beamforming, for example to manage interference 
between users. To multiplex several data streams and perform 
more accurate beamforming, hybrid precoding was proposed 
0 0 where the processing is divided between the analog 
and digital domains. In 0 0 the sparse nature of the 
mmWave channels was exploited to develop low-complexity 
hybrid precoding algorithms for single-user channels that 
support a limited number of streams 0. The digital precoding 
layer of hybrid precoding gives more freedom in designing the 
precoders, compared to the analog-only solution, and seems 
like a promising framework on which to build multi-user 
MIMO mmWave precoders. 

In this paper, we propose a low-complexity yet efficient 
hybrid analog/digital precoding algorithm for downlink multi¬ 
user mmWave systems. The proposed algorithm depends on 
the known (but arbitrary) array geometry and incurs a low 
training and feedback overhead. Our model assumes that the 
MS’s employ analog-only combining while the BS performs 
hybrid analog/digital precoding. We provide a lower bound on 
the achievable rate using the proposed algorithm for the case of 
single-path channels, which is relevant for mmWave systems, 
and provides useful insights for more general settings. The 
proposed algorithm and performance bounds are also evaluated 
by simulations and compared with analog-only beamforming 
solutions. The results indicate that the proposed algorithm 
performs very well in mmWave systems thanks to the sparse 
nature of the channel and the large number of antennas used 
by the BS and MS’s. An extended journal version of this 
work was submitted recently which provides an analysis 
of the proposed algorithm in more general channel settings, 
and characterizes the rate-loss when a limited feedback exists 
between the BS and MS’s. 

We use the following notation: A is a matrix, a is a vector. 



Fig. 1. A multi-user mmWave downlink system model, in which a BS uses 
hybrid analog/digital precoding and a large antenna array to serve U MSs. 
Each MS employs analog-only combining to receive the signal. 
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a is a scalar, and ^ is a set. ||A||i? is the Frobenius norm of 
A, whereas A^, A*, A“^, are its transpose, Hermitian, and 
inverse, respectively. E [•] denotes expectation. 

II. System Model 

Consider the mmWave system shown in Fig. A BS with 
Nbs antennas and N^f RF chains is assumed to communicate 
with U MS’s. Each MS is equipped with Nms antennas as 
depicted in Fig. We focus on the multi-user beamforming 
case in which the BS communicates with every MS via only 
one stream. Therefore, the total number of streams is As = U. 
Further, we assume that the maximum number of users that 
can be simultaneously served by the BS equals the number of 
BS RF chains, i.e., U < Arf. This is motivated by the spatial 
multiplexing gain of the described multi-user hybrid precoding 
system, which is limited by mm{NBF,U) for Nbs > ^rf- 
For simplicity, we will also assume that the BS will use U out 
of the Nbf available RF chains to serve the U users. 

On the downlink, the BS applies slUxU baseband precoder 
Fbb = •••5 followed by an Abs x U RF pre¬ 
coder, Frf = •••, • The sampled transmitted 

signal is therefore 

x = FrfFbbs, (1) 

where s = [51,82,..., 5^7]^ is the f/ x 1 vector of transmitted 
symbols, such that E [ss*] = and P is the average total 
transmitted power. Since Frf is implemented using quantized 
analog phase shifters, [Frf]„j^„ = where 

is a quantized angle, and the factor of is for power 

normalization. 

For simplicity, we adopt a narrowband block-fading channel 
model as in 0, 0, 0 in which the uth MS observes the 
received signal as 

u 

r„ = H„ y] FRpf^^Sn + n„, (2) 

n=l 

where is the Ams x Ars matrix that represents the 
mmWave channel between the BS and the uth MS, and 
Du ^ cr^I) is a Gaussian noise vector. 

At the uth MS, the RF combiner is used to process the 
received signal to produce the scalar 

u 

=w*H„yyFRFf®®s„+w*n„, (3) 

n=l 


Fig. 2. A BS with hybrid analog/digital architecture communicating with 
the uth MS that employs analog-only combining. 


where v^^^ has similar constraints as the RF precoders, i.e., 
the constant modulus and quantized angles constraints. We 
assume that only analog (RF) beamforming is used at the MS’s 
as they will likely need cheaper hardware with lower power 
consumption. 

MmWave channels are expected to have limited scatter¬ 
ing 0. To incorporate this effect, we assume a geometric 
channel model with Lu scatterers for the channel of user u. 
Each scatterer is assumed to contribute a single propagation 
path between the BS and MS 0- The channel model can be 
transformed into the virtual channel model, which simplifies 
the generalization for larger angle spreads by incorporating 
spatial spreading functions (Tj . Under this model, the channel 
can be expressed as 

H„ = ^ ^ ^4^ 

^ ^ e=i 

where auy is the complex gain of the path, including 
the path-loss, with E [|(Tn,^P] = The variables 0 u,£, and 
^ [0? 27 r] are the path’s angles of arrival and departures 
(AoA/AoD) respectively. Einally, aes ( 0 u/) and aMS 
are the antenna array response vectors of the BS and uth 
MS respectively. The BS and each MS are assumed to know 
the geometry of their antenna arrays. While the algorithms 
and results developed in the paper can be applied to arbitrary 
antenna arrays, we use uniform planar arrays (UPAs), with 
both azimuth and elevation components of each angle, in the 
simulations of Section [Vll 


III. Problem Eormulation 


Given the received signal at the uth MS in ([^, which is then 
processed using the RE combiner v^^^, the achievable rate of 
user u is 


Ru = log 2 1 + 


P|w*H„FRFfri 


FE„^JwJlH„FRFfBB|^ + <T2 


( 5 ) 


The sum-rate of the system is then i?sum = 

Our main objective is to efficiently design the analog (RE) 
and digital (baseband) precoders at the BS and the analog 
combiners at the MS’s to maximize the sum-rate of the system. 










































Due to the constraints on the RF hardware, such as the 
availability of only quantized angles for the RF phase shifters, 
the analog beamforming/combining vectors can take only 
certain values. Hence, these vectors need to be selected from 
finite-size codebooks. There are different models for the RF 
beamforming codebooks, two possible examples are 

1) General quantized beamforming codebooks Here, the 
codebooks are usually designed for rich channels and, 
therefore, attempt a uniform quantization of the space of 
beamforming vectors. These codebooks were commonly 
used in traditional MIMO systems. 

2) Beamsteering codebooks The beamforming vectors, 

here, are spatial matched filters for the single-path chan¬ 
nels. As a result, they have the same form of the array 
response vector and can be parameterized by a simple 
angle. Let represent the RF beamforming codebook 
with cardinality \J^\ = A^q. Then, T consists of the 
vectors aes variable /cq taking the values 

0,1, 2, and Aq — 1. The RF combining vectors codebook 
W can be similarly defined. 

Motivated by the good performance of single-user hybrid 
precoding algorithms 0 0 which rely on RF beamsteering 
vectors, and by the relatively small size of these codebooks 
which depend on a single-parameter quantization, we will 
adopt the beamsteering codebooks for the analog beamform¬ 
ing vectors. While the problem formulation and proposed 
algorithm in this paper are general for any codebook, the 
performance evaluation of the proposed algorithm done in 
Section |V| depends on the selected codebook. 

If the system sum-rate is adopted as a performance metric, 
the precoding design problem is then to find F^p, 
and that solve 


{Frf,{C®®}Li >{<}«=! } 


arg max y: l0g2 1 + 


P|w:H„FRFfBBr 


U=1 


FE„^Jwj;H„FRFfBB|^ + a2 


s.t. [Frf].^^ G a, rt = 1, 2,..., U, 

IIFkp [ff^f,^^...,f^^] 111 = A 


( 6 ) 


The problem in ([^ is a mixed integer programming prob¬ 
lem. Its solution requires a search over the entire x 
space of all possible Frf and combinations. Further, 

the digital precoder Frb needs to be jointly designed with 
the analog beamforming/combining vectors. In practice, this 
may require the feedback of the channel matrices = 

1,2, or the effective channels, w^H^^Frf. Therefore, 

the solution of (|^ requires large training and feedback over¬ 
head. Moreover, the optimal digital linear precoder is not 
known in general even without the RF constraints, and only 
iterative solutions exist ©• Hence, the direct solution of this 
sum-rate maximization problem is not practical. 


Algorithm 1 Multi-user Hybrid Precoders Algorithm 


Input: T and W, BS and MS RF beamforming codebooks 
First stage: Single-user RF beamforming/combining design 
For each MS u,u = 1, 2,..., U 

The BS and MS u select v* and g* that solve 
{§«><} = argmax||g;H„v„|| 


MS u sets vju = 
BS sets Frf = [v|,v: 


Vg^ew 


•• 


Second stage: Multi-user digital precoding design 
For each MS u,u = 1, 2,..., U 

MS u estimates its effective channel = w^H^^Frf 
MS u feeds back its effective channel h,, to the BS 


BS designs Fbb = H 



r—T —T 1 

(hH j , H = 



1 T 




l|FHFfBB|| 


,u= 1,2,...,P 


Similar problems to ([^ have been studied before in litera- 
ture, but with baseband (not hybrid) precoding and combining 
0-(n) The application of these algorithms in mmWave 
systems, however, is generally difficult due to (i) the large 
feedback overhead associated with the large antenna arrays 
at both the BS and MS’s, (ii) the need for the RF beam¬ 
forming/combining vectors to be taken from quantized code¬ 
books which requires a search over both the BS and MS’s 
beamforming vectors and limits the freedom in designing 
the combining vectors, and (iii) the convergence of iterative 
precoding/combining algorithms like 0 has not yet been 
studied for hybrid precoders. 

Given the practical difficulties associated with applying the 
precoding/combining algorithms in mmWave systems, we pro¬ 
pose a new mm Wave-suitable multi-user MIMO beamforming 
algorithm in Section W Our proposed algorithm is developed 
to achieve a good performance compared with the solution of 
0, while requiring low training and feedback overhead. 


IV. Two-stage Multi-user Hybrid Precoding 

The additional challenge in solving (7), beyond the usual 
coupling between precoders and combiners 0-|[^, is the 
splitting of the precoding operation into two different domains, 
each with different constraints. To overcome that, we propose 
Algorithm that divides the calculation of the precoders into 
two stages. In the first stage, the BS RF precoder and the MS 
RF combiners are jointly designed to maximize the desired 
signal power of each user, neglecting the resulting interference 
among users. In the second stage, the BS digital precoder is 
designed to manage the multi-user interference. The operation 
of Algorithm can be summarized as follows. 

In the first stage: The BS and each MS u design the 
RF beamforming and combining vectors, f^^ and to 
maximize the desired signal power for user u, and neglecting 
the other users’ interference. As this is the typical single-user 
RF beamforming design problem, efficient beam training al¬ 
gorithms developed for single-user systems that do not require 
explicit channel estimation and have a low training overhead. 











can be used to design the RF beamforming/combining vectors. 
For example, 0, © developed multi-resolution codebooks 
to adaptively refine the beamforming/combining vectors, and 
hence avoid the exhaustive search complexity. 

In the second stage: The BS trains the effective channels, 
hu = = 1,2,with the MS’s. Note that 

the dimension of each effective channel vector is U xl which 
is much less than the original channel matrix. Then, each MS 
u feeds its effective channel back to the BS, which designs 
its zero-forcing digital precoder based on these effective 
channels. Thanks to the narrow beamforming and the sparse 
mmWave channels, the effective channels are expected to be 
well-conditioned which makes adopting a multi-user digital 
beamforming strategy like zero-forcing capable of achieving 
near-optimal performance as will be shown in Section [V] 

V. Achievable Rate with Single-Path Channels 

The analysis of hybrid precoding is non-trivial due to the 
coupling between analog and digital precoders. Therefore, we 
will study the performance of the proposed algorithm in the 
case of single-path channels. This case is of special interest 
as mmWave channels are likely to be sparse, i.e., only a 
few paths exist Further, the analysis of this special case 
will give useful insights into the performance of the proposed 
algorithms in more general settings. 

Next, we will characterize a lower bound on the achievable 
rate by each MS when Algorithm[^is used to design the hybrid 
precoders at the BS and RF combiners at the MS’s. Consider 
the BS and MS’s with the system and channels described in 
Section |I^ with the following assumptions: 

Assumption 1: All channels are single-path, i.e., = 1, 

u = For ease of exposition, we will omit the 

subscript I in the definition of the channel parameters in Q. 

Assumption 2: The RF beamforming and combining vec¬ 
tors f^^ and v^n,^ = l,2,...,f/, u = l,2,...,f/ are beam¬ 
steering vectors with continuous angles. 

Assumption 3: The BS perfectly knows the effective chan¬ 
nels hu^ 14 = 1, 2,..., U. 

In the first stage of Algorithm the BS and each MS u 
find V* and g* that solve 

{g«>v*} = argmax||g;H„v„||. ( 7 ) 

vg„ew 

As the channel has only one path, and given the 
continuous beamsteering capability assumption, the optimal 
RF precoding and combining vectors will be g* = aMs(^u), 
and V* = aBs(0u)- Consequently, the MS sets = aMs(^u) 
and the BS takes f^^ = aBs(0w)- If we let the A^bs x U matrix 
Abs gather the BS array response vectors associated with the 
U AoDs, i.e., Abs = [^bs (0i) , ^bs (^ 2 ), •••, ^bs we 

can then write the BS RF beamforming matrix, including the 
beamforming vectors of the U users as Frf = Abs- 

The effective channel for user u after designing the RF 
precoders and combiners is 

= WuUu'Prf 

= VNbsJ^MSO^u^bS i^u) Frf- 


_ _ ^ ^ ^ 

Now, defining H = [h^, h 2 ,h^]^, and given the design 
of Frf, we can write the effective channel matrix H as 

H = DA*sAbs, (9) 

where D is a f/ x f/ diagonal matrix, [D]^ ^ = V^bs^msc^u- 
Based on this effective channel, the BS zero-forcing digital 
precoder is defined as 

Fbb = H* ( 10 ) 

where A is a diagonal matrix with the diagonal ele¬ 
ments adjusted to satisfy the precoding power constraints 
||FRFf^^||^ = 1,14 = 1,2,...,[/. The diagonal elements of 
A are then equal to [See Appendix A in |[^ for a derivation] 




NbsNms 


(Ars^BS 


,u = l,2,...,f/. 


( 11 ) 


Note that this A is different than the traditional digital zero¬ 
forcing precoder due to the different power constraints in the 
hybrid analog/digital architecture. 

The achievable rate for user 14 is then 
SNR 



U 


h f 


BB 


SNRiVBsiVMs|«ur 
U (A^sAbs)”^ 


( 12 ) 


To bound this rate, the following lemma, which character¬ 
izes a useful property of the matrix AggAes. can be used. 


Lemma 4: Let Ars = [^bs (0i) , ^bs (^ 2 ), ^bs (0i/)], 
with the angles (j)u^u = l,2,...,f/ taking continuous values 
in [0,27r], then the matrix P = AggARs is positive definite 
almost surely. 

Proof: Let the matrix P = AbsAbs, then for any non¬ 
zero complex vector z G it follows that z*Pz = 

IIAbszIII > 0. Hence, the matrix P is positive semi-definite. 
Further, if the vectors aRs (0i), aRS (^ 2 ), ...,aBs(0i/) are 
linearly independent, then for any non-zero complex vector z, 
Absz / 0, and the matrix P is positive definite. To show that, 
consider any two vectors aRs (fu ), aRs (0n)- These vectors 
are linearly dependent if and only if fu = fn- As the 
probability of this event equals zero when the AoDs fu and 
fn are selected independently from a continuous distribution, 
the matrix P is positive definite with probability one. ■ 
Now, using the Kantorovich inequality fT^ , we can bound 
the diagonal entries of the matrix (AggARs)”^ using the 
following lemma from GD 

Lemma 5: For any n x n Hermitian and positive definite 
matrix P with the ordered eigenvalues satisfying 0 < Amin < 


A 2 < ... < Amax, the element (P)^{^ ,i4 = 1,2, ...,n satisfies 


(p)“^ < ^ 

V /w,u — 


Amax (^) I ^min (P) 

4[P]n,u V ^min (P) Amax (P) 


(13) 


We also note that for the matrix AggARs, we have 
(Abs-^bs)„,„ = h Amin (Aj^gABs) = (Abs), and 















Fig. 3. Achievable rates per-user using the hybrid precoding and beamsteer¬ 
ing algorithms with perfect channel knowledge. 



Fig. 4. Achievable rates per-user using the hybrid precoding and beamsteer¬ 
ing algorithms for different values of angle spread. 


Amax(A^sABs) = cr^ax Abs), where CTi„ax(ABs) and 
c^min(ABs) dfe the maximum and minimum singular values, 
respectively. Using this note and Lemma and defining 
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we can bound the achievable rate of user u in as 


(14) 


Ru > log 2 (l + ^AbsVms \au\^G , (15) 


In addition to characterizing a lower bound on the rates 
achieved by the proposed hybrid analog/digital precoding 
algorithm, the bound in ( p?] ) separates the dependence on the 
channel gains au, and the AoDs = which can 

be used to claim the optimality of the proposed algorithm in 
some cases and to give useful insights into the gain of the 
proposed algorithm over analog-only beamsteering solutions. 
This is illustrated in the following results. 

Proposition 6: Let Ru = log 2 (^1 + ^AbsAms 
denote the single-user rate, and let denote the rate 

achieved by user u when the BS employs analog-only beam¬ 
steering designed based on the first stage of Algorithm 


log2 (l 


^JVBsiVMs|a„|" 




)■ 


with 


Pu,n = Sbs (</>«) aBs (</>n)- When Algorithm is used to 
design the hybrid precoders and RF combiners described in 
Section |n| and given the Assumptions 1-3, the achievable rate 
by any user u satisfies 


1) E 


Ru — R 


2) limATss ^oo 

3) limiVMs ^OO 


uj < A(Vbs,(/), 

Ru = Ru with probability one, 
E [Ru - i?®®] = oo, 


where K (Abs, U) is a constant whose value depends only on 
Abs and U. 

Proof: Please refer to Appendices B and C in j^. ■ 


Proposition indicates that the average achievable rate 
of any user u using the proposed low-complexity precod¬ 
ing/combining algorithm grows with the same slope of that 
of the single-user rate at high SNR, and stays within a 
constant gap from it. This gap, A(ABs,f^), depends only 
on the number of users and the number of BS antennas. As 
the number of BS antennas increases, the gap between the 
achievable rate using Algorithm and the single-user rate 
decreases, and approaches zero at infinite antenna numbers. 
One important note is that this gap does not depend on the 
number of MS antennas, which is contrary to the analog-only 
beamsteering, given by the first stage only of Algorithm 
This leads to the third part of the proposition. This implies 
that multi-user interference management is still important at 
mmWave systems even with large numbers of antennas at the 
BS and MS’s, and perfect alignment. Note also that this is not 
the case when the number of BS antennas goes to infinity as it 
can be easily shown that the performance of RF beamsteering 
alone becomes optimal in this case. 

VI. Simulation Results 

In this section, we evaluate the performance of the proposed 
hybrid analog/digital precoding algorithm and the presented 
bounds using numerical simulations. 

First, we compare the achievable rates without quantization 
loss in Fig. where we consider the system model in 
Section \n\ with a BS employing an 8 x 8 UFA with 4 MS’s, 
each having a 4 x 4 UFA. The channels are single-path and 
Rayleigh distributed, the azimuth and elevation AoAs/AoDs 
are assumed to be uniformly distributed in [0, 27r] and [— f, f ], 
respectively. The rate achieved by the proposed algorithm is 
compared with single-user and beamsteering rates. The figure 
shows that the performance of hybrid precoding is very close 
to the single-user rate thanks to canceling the residual multi¬ 
user interference. The gain over beamsteering increases with 
SNR as the beamsteering rate starts to be interference limited. 






































The tightness of the derived lower bound is also shown. 

To evaluate the performance of the proposed hybrid pre¬ 
coding algorithm in more general channel settings and system 
imperfections, we consider a 3-cluster channel model in Fig.|^ 
Each cluster has 6 rays with Laplacian distributed AoAs/AoDs 
0 - The phase shifters at the BS and MS’s are assumed to 
be quantized with 6 bits and 4 bits, respectively. The figure 
shows that when angle-spread and system imperfections exist, 
the performance of hybrid precoding is still within a small gap 
from the single-user rate. Further, it provides a good gain over 
analog-only solutions. 

Finally, Fig. evaluates the performance of the proposed 
algorithm in a mm Wave cellular setup including inter-cell 
interference, which is not explicitly incorporated into our de¬ 
signs. In this setup, BS’s and MS’s are assumed to be spatially 
distributed according to a Poisson point process with MS’s 
densities 30 times the BS densities. The channels between the 
BS’s and MS’s are single-path and each link is determined 
to be line-of-sight or non-line-of-sight based on the blockage 
model in j^. Each MS is associated to the BS with less path- 
loss and the BS randomly selects n = 2, ..,5 users of those 
associated to it to be simultaneously served. BS’s are assumed 
to have 8x8 UPAs and MS’s are equipped with 4x4 UPAs. 
All UPA’s are vertical, elevation angles are assumed to be 
fixed at 7r/2, and azimuth angles are uniformly distributed in 
[0, 27r]. Fig. shows the per-user coverage probability defined 
as 7^ (Ru > r/), where r] is an arbitrary threshold. This figure 
illustrates that hybrid precoding has a reasonable coverage gain 
over analog-only beamsteering, especially when large numbers 
of users are simultaneously served, thanks to the interference 
management capability of hybrid precoding. 
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VIII. Conclusions 

In this paper, we proposed a low-complexity hybrid 
analog/digital precoding algorithm for downlink multi-user 
mmWave systems inspired by the sparse nature of the channel 
and the large number of deployed antennas. For the single-path 
channels, we characterized the achievable rates of the proposed 
algorithm and illustrated its asymptotic optimality. Compared 
with analog-only beamforming solutions, we showed that 
higher sum-rates can be achieved using hybrid precoding. The 
results indicate that interference management in multi-user 
mm Wave systems is still important even when the number 
of antennas is very large. Simulation results verify these 
conclusions and show that they can be extended to more 
general multi-path mmWave channel settings. 
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